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ABSTRACT 

The radiation from the central regions of active galactic nuclei, including that from the 
accretion disk surrounding the black hole, is likely to peak in the extreme ultraviolet ~ 13 — 
100 cV. However, due to Galactic absorption, we are limited to constraining the physical 
properties - black hole mass and accretion rate - from what observations we have below ^ 
10 eV or above ~ 100 cV. In this paper we predict the thermal and ionization states of warm 
absorbers as a function of the shape of the unobservable continuum. In particular we model 
an accretion disk at kXm ^ 10 cV and a soft excess at kTse ^ 150 eV. The warm absorber, 
which is the highly ionized gas along the line of sight to the continuum, shows signatures 
in the ~ 0.3 — 2kcV energy range consisting of numerous absorption lines and edges of 
various ions, some of the prominent ones being H- and He-like oxygen, neon, magnesium and 
silicon. We find that the properties of the warm absorber are significantly influenced by the 
changes in the temperature of the accretion disk, as well as by the strength of the soft excess, as 
they affect the optical depth particularly for iron and oxygen. These trends may help develop 
a method of characterising the shape of the unobservable continuum and the occurrence of 
warm absorbers. 
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1 INTRODUCTION 

Accretion of matter onto a supermassive black hole and the 
conversion of the gravitational energy into radiation, via an 
accretion disk, is the main source of energy output in ac- 
tive galactic nuclei (AGN). The radiation from AGN is likely 
to peak in the ex tr eme ultraviolet (EUV) ~ 10 — 100 eV 
jLvnden-Belll 1 19691: [Pringle. Rees & Pacholczvkl 1 19731 : IShieldj 
1 19781 : IShakura & Sunvaevlll973l ; lNetzeill 19851) . However. Galactic 
absorption of EUV light results in a 'blind spot' in the energy range 
13 — 100 eV (~ 912 — 50A, marked as "unobserved" in Figure[TJ 
introducing an uncertainty in the shape of the spectral energy distri- 
bution (SED) of AGN, just where it is expected to peak. Attempts 
have been made to reconstruct the SED fro m observations avail- 
able below ~ 10 eV and above ~ 100 eV jSiemiginowska et al. I 
Il995l :[ Sobolewska et al. ,l2004al lhl and references therein). This un- 
observed part of the SED is most effective at ionizing gas (which is 
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the very reason it also gets absorbed by the Galaxy). Hence the con- 
tinuum will influence the intervening line of sight gas, local to the 
AGN, and the resultant absorption and emission line strengths of 
the various ions seen in far ultraviolet (EUV) and the Balm er con- 
tinuum can be used to constrain the SED, as discussed bv iNetzeJ 
( Il987l) . 

At the high energy end of the 'blind spot', i.e. in the soft X-ray 
(0.1 - 10 keV), observations of most AGN reveal that, for a major- 
ity, the SED from 2-10 keV is well approximated by a power-law 
f{i') ~ with spectral indices ~ 0.8. The join between the 
EUV (2500A) and the soft X-ray (2 keV), apx, is steep er, with 
typical values of ~ 1.2 for Seyfert 1 galaxies l lNetzeill993h . Hence 
extending the X-ray power-law to the ultrav iolet (UV) cannot ex- 
plain the observed flux (but see lLaorl l ll997l) ). An ionizing contin- 
uum with two power-law components, a steep one for the join be- 
tween EUV and soft X-ray, and a flatter one for the soft X-ray SED, 
can represent the overall shape of the AGN spectra in the energy 
range 3 eV — 2 keV. 

However, such an SED is not a physical model and does 
not resemble certain specific features expected in this energy 
range. In particular, the proposed accretion disk, whose in- 
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nemost stable orbit should have a temperature of ~ 13 eV, 
for a black hol e of mass 10^ and an accretion rate of 
rn/ir iEdd ~ 0.1 jShakura & Sunvaevlll973l ; lFrank. King & Raind 
is not included. Another spectral component which 
would remain unaccounted for by a simple two power-law 
continuum is the soft excess which is seen in most AGN below 
2keV. The soft exces s coul d be a Comptonised disk compo- 



nent I Czemv & ElvisI 



---- - - | l987h. an additional quasi-blackbod' 

dCmmmv et al. I 120061: iKorista. Ferland & BaldwiiJ 
20051) 



3; 



20051) or due to relativistic broadened 



IRoss &Fabianril993r 

absorption features jGierlinski & Pond |2004|) . Such spectral 
components are likely to have important effects on the nature of 
the absorbing gas along the line of sight towards the central engine 
of the AGN. 

In the 0.3 - 1.5 keV soft X-ray spectra of about half of 
all Seyfert 1 galax ies jNandra & Pounds! 1 1994 IRevnoldsl [19971 : 
lGeorgeetal.il 19981) and quasars jPiconcelli et al. Il2005l) ~one can 
find signatur es of pho t oioniz ed gas called "warm absorbers" (here- 
after WAs, IHalpeml ( Il984h ). These signatures are, often, ab- 
sorption lines and edges from highly ioni zed species, such as 
OVII, OVIII, FeXVII, NeX, CV and CVI jColUnge et al. Il200ll : 



Kaastra et al. "2002'; 'Kinkhabwala et al. "2002':'Blustin et al. l2003l : 
kron gold et al . 2003; Netzer et al. , 2003 ; Tur ner e t al. j |2004|) . 



The typical column density observed for the gas is A^h ~ 
■|^q22±i jjjg ionization parameter ^ of the WA spans a range 

~ 10 — 1000 ergcms"^ corresponding to gas temperatures of 
10** — 10^'^ K. Some authors claim that the WA is made up 
of discrete thermodynamic phas es in near pressure equilibrium 
jAndrade- Velazquez et al. 2010l and references therein), w hile 



others jOgle et al. l2004l : ISteenbrugge et al. l2005l ; lBehaj2009l) be 

lieve that the gas has a continuous distribution of temperature and 
pressure. 

In this temperature range heating and cooling processes are 
respectively dominated by photoionization and recombination of 
higher ionization states of heavier elements. These processes re- 
sult in local region s of thermal stability in this otherwise unstable 
temperature rane e jKrolik et al. Ill98ll: iGehrels & Williamsll 19931 : 
lHessetal.il 19971) . IChakravortv et al. I ilOO'^ showed that the na- 
ture of the WA is strongly influenced by the chemical composition 
of the absorbing gas. It was shown that the abundance of iron and 
oxygen, which have important atomic transitions in the sub-keV 
energy range are particularly important. 

To obtain the equilibrium conditions for a gas with an assumed 
set physical conditions, we have used the publicly available pho- 
toionization code CLOUD"'i0version C07.02 (hereafter C07), see 
Ferland (1998). CLOUDY calculates the ionization equilibrium 
conditions by solving the energy and charge conservation equa- 
tions under the assumption that all the atomic processes have had 
time to reach a steady state. Any stable photoionised gas will lie 
on the thermal equilibrium curve (hereafter referred to as the sta- 
bility curve) where heating balances cooling. The stability curve 
of temperature (T) against pressure ((,/T), where ^ is the ioniza- 
tion parameter (see below for definition), is often used to study 
the nature of the WA. Gas lying off the stability curve will heat 
or cool until reaching the curve. If the curve has kinks that pro- 
duce multiple stable values at fixed £,/T, then the WA can have 
multiple temperature phases in pressure equilibrium. The shape 
of the stability curve depends on the SEP of the ionizing contin- 
uum and the chemical abundance of the gas CRevnolds & Fabiaa 



1993: iKrolik & Krissll200ll ; iKomossa & MeerschweinchenI 1200(3: 



Komossa & Mathuj|200ll : IChakravortv et al. Il2008ll2009l) . 

In this paper we investigate the behaviour of the stability curve 
as a function of the disk blackbody and soft excess both to under- 
stand the role of these spectral components on the nature of the WA, 
and to examine whether the observable parameters of the WA can 
diagnose the shape of the AGN continuum in the EUV 'blind spot'. 
In Section |2| we describe the different models of the AGN contin- 
uum that we are interested in and use these SEPs to study their 
effect on the stability curves in Section|3] We investigate the causes 
of the variations in the stability curves in Section|4|by studying the 
heating and ionization fractions of the elements and ions which are 
responsible for determining the nature of the WA. Section|5]exam- 
ines the relevance of a few other physical parameters, like the X-ray 
slope of the ionizing conti nuum and the abundance of the absorb- 
ing gas, (also discussed in IChakravortv et al. II2009I) in the context 
of the present study. The multi-phase nature of the WA is discussed 
in Section|6|and we conclude our results in Section|7| 



2 THE AGN CONTINUUM 

The 2-10 keV X-ray spectra of most of the observed active galac- 
tic nuclei (AGN) can be modeled satisfactorily using a power-law 
/{v) ~ u^" where the spe ctral index a for mos t AGN lie in 
the range 0.7 < a < 0.9 jWilkes & ElvisI Il987l : ICrupe et al. I 
I2OO6I : iLopez et al. II2OO6I) . The optical depth r = 1 for photons 
at 0.2 keV with a typical Galactic hydrogen column density of 
Nu = 3 X 10^° cm~^. Extinction is a steeply rising function of 
column density, e. g. r = 2 at 0.2 keV for iVn = 5 x 10^" cm"^. 
Thus, we lack observational constraints on the evolution of the X- 
ray power-law at lower energies. However, for most AGN, the UV 
flux is found to be much higher than what a simple extr apolation of 
the X-ray power-law to lower energies woul d predict jElvis et alTI 
ll994IZheng et al. II 19971 : IShang et al. II2005I) , but see lUaoJ l ll99'7h 
for exceptions. The slope aox, defined as 



aox = -0.384 log 



fi2keV) 



(1) 



_/(4.7eV = 2500A) 

bv lTananbaum et aTl ( 1 19791) . is conventionally used to parametrize 
the nominal power-law between th e UV and soft X-ray bands. The 
observed range is 1 < aox < 2 l lStalin et al. II2009I : iGreen et al. I 



12003, ™d references therein). 

The overall shape of the AGN continuum in the energy range 
of ~ 4.7 eV(2500A) — 10 keV can, thus, be represented as a sum 
of two power-law components 



for E(= hu) > 4.7 eV, (2) 



where as{> 0) is the spectral index of a steep soft component and 
Al is the relative normalization factor to realiz e the desired values 
of ao x (dotted-and-dashed curve in FigurefTI). IChakravortv et 
i2009*) had used such an SEP with Er^ax = hiymax = 200 keV, 
and this value is maintained throughout this paper, as well. Results 
from the Swift/BAT (Burst Alert Telescope) hard X-ray sky survey 
jTueller et ani2008h show that, of the brightest few dozen AGN, 
for which Em ax can be determined, the cut-off spans the range 50 
to 450 keV. In lChakravortv et alTI ^20091) we have shown that in the 
range 50 keV < Emax < 400 keV, the variation in Emax does not 
have any significant effect on the nature of the warm absorber. For 
energies lower than 4.7 eV IChakravortv et al. I ( |2009|) followed the 
cut-off scheme 



^ URL: |http://www.nublado.org/] 



for 2.8 < E(= hv) < 4.7 eV 
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Figure 1. Comparison between various ionizing continua : the broken 
power-law (Equation |2] dotted-and-dashed line) and more realistic SEDs 
having disk blackbody and soft excess components (Equations|4]and|5). For 
all the curves a. = 0.8 and aox = 1-2. Using Equation |4] we draw the 
solid curve with rj = 0.0 and the dashed curve with rj = 1.0, while the 
dotted curve is drawn using Equation [s] with rj' = 1.0. The temperature of 
the innermost ring of the accretion dislc is Ti„ = 10 eV and temperature 
corresponding to the soft excess is Tse = Tc = 150 eV. See text for further 
details. 



such, the resultant shape of the "Big Blue Bump" peaks between 
~ 10 — 100 eV for the range of black hole masses and accretion 
rates typical of Seyfert galaxies. 

Studies of high redshift quasars show that their UV-EUV 
SEDs are consistent with the standard th eor y of accretion disk s 
teechtold et all 1 1994 IPage et al. I l2004allbl (shang et al. 1 12005[) . 
However, typ ically, the quasars have higher mass (10^ < Mbh 
5 X 10"'M,7, [siiang et al. II 20051: iBechtold et al. II 19941) and higher 
accretion rates teechtold etal. 1 1994h . Thus, for the same accretion 
rates, we can expect the "Big Blue Bump" in quasars to be peaking 
at lower energies than th at in Seyferts which w ould typically result 
in steeper aox ^ 1- 8 teechto ld et al. |[T994|) for the quasars, as 
compared to 1.2 jNetz er 1993) for the Seyferts. 

A standard model for the spectral co mponent from the ac- 
cretion disk is availabl e as disk bla ckbody I Mitsuda et al. 1 1 19841 : 
iMakishima et al. II 19871) in XSPEC0 ( lArnaujll996l) . We have used 
version 11.3 of XSPEC to generate the disk blackbody spectral 
component fdbb{i') to be used in the construction of the realistic 
AGN continuum (the third term in Equations |4| and O. The solid 
black line in Figure [T] shows a SED which includes a "Big Blue 
Bump" with Tin = 10 eV. In this paper we explore "Big Blue 
Bumps" in the range 10 < r,„ < 30 eV corresponding to lO* "*^ < 
A/sH < Vf''^^ Mq for an accretion rate of m/mEdd = 0.1. This 
range of accretion disks parameters pertain to Seyfert galaxies or 
low mass and low accretion rate quasars. We do not explore the pa- 
rameter range for higher mass quasars because for them Tin would 
move to lower energies and such disk blackbodies are less likely to 
influence the nature of the WA. 



~ ° for 0.12 < E(= hi/) < 2.8eV 

~ u^-^ for E(= hi^) < 0.12 eV (3) 

This s cheme is similar to that described by iMathews & Ferlan3 
( ll987l) . However, a SED described by Equations|2|and|3]may some- 
times be an inadequate description for the AGN ionizing continuum 
because it fails to represent the signature of the accretion disk at 
~ 10 eV or the soft excess component. 



2.1 Disk blackbody 

Multiwavelength observations suggest that the AGN continua peak 



the Quasar luminositv jNeugebauer et al. 1 1 19791: IShieldsl 


197S: 


Malkan & Sargent 19821: 


Elvis et al. 1 


19861: ICzernv& El vid 


1987: 


Mathews & Ferlandl 1987 


: Laor et al. 


1990l: Siemiginowska et al. 


19951: IZheng et al. II 19971: ISobolewska et al. Il2004dlhl). This soec- 



tral component, often referred to as the "Big Blue Bump", is con- 
sidered to be th e signature of the presence of an accretion disk 
as dis c ussed by 



( ll973l) : IShieldsi 



s signature or tne p r esence ot an accretion disk 
Lvnden-Belll ([196 91): IPringle. Rees & Pacholczvkl 
1978l) : IShakura & Sunvaevj h973l) . 



According to the st andard theory of accretion disks by 
IShakura & Sunvaevl l ll973h . the emission from the disk can be mod- 
eled as a sum of local blackbody radiations emitted from annuli of 
the disk at different radii. The temperature of all the subsequent 
rings can be estimated from the theory if the temperature Ti„ of 
the innermost ring of the accretion disk is known, which in turn is 
related to the one-fourth power of the mass of the central black 
hole and to the fourth power of its accretion rate. Thus for the 
same accretion rate, Ti„ changes by only a factor of 3 while the 
black hole mass spans a range of about two order of magnitude. As 



2.2 Soft excess 

X-ray observations of AGN with ROSAT and XMM-Newton 
often show that if the 1-10 keV power-law is extended to 
lower energies to fit the observed spectra of type 1 AGN, 
some unacc ounted for excess intensity is usually seen at 
E < IkeV (lElvis. Wilkes &Tananbauidll985l : lBrinkmaniJll9"93: 
iBuehler et al. I Il995l : IPounds & ReevesI l2002h . This excess has 
come to be known as the soft excess component. A blackbody with 
temperature Tse ~ 100 - 200 eV (i.e. pea king at ~ 282 - 564 eV) 



is often a good fit to the soft excess (iMatsumoto et al. 
IPiconcelli et aTl l2005l : IPorguet et al"! |2004|: I Vignali et al. 



2004 



2004 



and references therein). The same authors show that the ratio of 
the soft excess luminosity to power-law luminosity, usually be- 
twee n 0.1 - 10 keV varies fr om object to object fro m 0.04 in Mkn 
304 (IPiconcelli et al. I [200^ to > 1.0 in Ark 564 dVignali et alTI 
I2OO4I) . The simple 'sum of blackbodies' model for the accretions 
disk is a satisfactory qualitative representation of the UV SED with 
a range of parameter values to cover the observed AGN proper- 
ties. However, all AGN disks, having supermassive black-holes at 
their centres, are too cold to reach soft X-ra ys at ~ 0.5 keV . More 
sophisticated modifications of thi s model dCzemv & Elvis|[l987l : 
iKorista. Ferland & Baldwin! 1 19971) or an additional spectral com- 
ponent is required to explain the soft excess. 



2.2.1 Blackbody soft excess 

A theoretical representati on of the soft excess as a bl ackbody is 
a simplified version of the IRoss &FabiaiJ ( fl99ll2005l) model for 
the soft excess component which owes its origin to the reflection 
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of the power-law component of the AGN spectrum from the ac- 
cretion disk. Thus, in this model the X-ray power-law photons are 
reprocessed, instead of the ones from the accretion disc. 

The general SED for the ionizing continuum including a 
blackbody soft excess component can be given as 



f{u) [{U- 



+ n- 



2nh 



exp{hi//kBTs, 

+ A2fdbb{l^,Tin) j I 



} 



(4) 



The first term in the above equation represents the X-ray power-law 
with spectral index a. The second term is the soft excess compo- 
nent, as a blackbody distribution of photons where Tse = 150 eV 
is the temperature of the blackbody and 77 determines the ratio of 
luminosity in the soft excess component to that in the power-law 
between 0. 1 - 10 keV. The third term in Equation |4] is the disk 
blackbody component which is parametrised by Tin, the temper- 
ature of the innermost ring of the accretion disk. A2 is the nor- 
malisation factor to attain the desired value of aox- Unless oth- 
erwise mentioned, throughout this paper, we have used q = 0.8 
and aox = 1.2, Emax = humax ~ 200 keV. For numerical con- 
venience we put a cut-off for /(f) at E = hv = 0.12 eV below 
which the flux drops off as u^'''. 

In Figure [T] both, the black solid SED and the blue dashed 
SED are drawn using Equation|4] For both SEDs, Tin = 10 eV so 
that the peak of the diskblackbody lies at ~ 30 eV. The solid line 
represents a SED which does not have any soft excess component, 
i.e.,?; — 0. On the other hand, the blue dashed SED with rj = 
1.0 has a moderately strong soft excess. See further discussions in 
Section l33TT] 



2.2.2 Comptonisation by ^ 150 eY plasma 

iKorista. Ferland & Baldwinl ( 1 19971) discuss the puzzle that many 
quasars which show soft ionizing continuum, would have an in- 
sufficient number of photons at the 54.4 eV He II edge to gen- 
erate the observed strength of He II emission seen in the same 
objects. To explain the strength of the He II line in Mrk 335, 
iKorista. Ferland & Baldwinl ( Il997h invoked a "double peaked" UV- 
EUV continuum where the second "bump" peaks at ~ 54 eV and 
contains energy comparable to the classical UV bump at lower en- 
ergies ('^ 10 eV). The 

A model, where the photons from the accretion disk are up- 
scattered by inverse Comptonization due to energetic electrons, 
does a better job of qualitatively satisfying the required strength of 
the UV-EUV SED for objects like Mrk 335. ICzernv & ElvisI jl987l) 
have shown that the models which account for the electron scatter- 
ing of the accretion disk photons are better than the simple 'sum of 
blackbodies' model. Similar thermal Comp tonizati on models hav e 
been worked o ut bvlUghtman & Zdziarskii (ll987h:ICoDPil ( Il992h ; 
iHaardtl ( ll993h : ICopDil ( Il999t) and lBeloborodovl ( Il999h among oth- 
ers. 

The thermal Compt onization mo del of'Lightman & Zdziarskii 

lll987 ) has b een used bv IZdziarski et al. ( 1996) and exte nded by 



Zvckil ( ll999l) and their model is included in XSPEC 12.5 l lAmaudI 
1996h as nthcomp. We have used nthcomp, where we have assumed 



that the seed photons coming from the accretion disk (modeled as 
disk blackbody) are reprocessed by the thermal plasma to generate 
sufficient photons at sub-keV to mimic the soft excess, usually ob- 
served in AGN. The high energy cut-off for the resulting soft excess 
feature is parametrised by the electron temperature Tc, whereas the 
low energy rollover is dependent on the effective temperature of 



the seed photons from the accretion disk, which in this case is 
parametrised by Tn- Between the low and high energy rollovers 
the shape of the spectrum is not necessarily a power law, but can be 
parametrised by an asymptotic power law index Tc which would 
physically be determined by the combination of electron scatter- 
ing optical d epth and electron tempera t ure, i.e. by the Com pton 
y-parameter jRvbicki & Lightmanlll986l) . lBeloborodovl ( 1 19991) had 
shown that there is a simple relation Fc « between the 

power-law index and the Compton y-parameter. Thus, in our model, 
Fc, used as an input, gives a measure of the extent of Compton re- 
processing; i.e. larger the value of Fc, lesser is the number of pho- 
tons reprocessed from the disk blackbody component to the high 
energy photons at ~ Tc. 

The SED with this alternative model for the soft excess fc{i^), 
generated by nthcomp, can be written as 



+ A2fdbb{l^,Tin) j e~ "ma:! , 



(5) 



where 77' is the ratio of luminosity of the soft excess component 
for E{— hv) > 0.3 keV to the luminosity of the power-law 
component in the energy range 0.1 — lOkeV. The orange dotted 
curve in Figure [T] is drawn using Equation [5] with Tn = 10 eV, 
Tc = 150 eV 77' = 1.0 and Fc = 1.7. Note that at £ = 
54.4 eV (i.e. log v = 16.12) the orange dotted SED has 1.4 times 
higher flux than the blue dashed SED due to a blackbody soft ex- 
cess and this factor grows to 13.1 at E = 100 eV. Thus SEDs like 
the orange dotted line are better suited to explain observations of 
objects like Mrk 335 or Mrk 478. See further discussions on this 

issu e in Section | 3.3.2l 

IChakravortv et al. I ( l2009h showed that the WA temperature 
range, 10*'^ < T < 10®'^ K, is strongly influenced by the chemi- 
cal composition of the absorbing gas, particularly by the abundance 
of iron and oxygen which have important atomic transitions in the 
sub-keV energy range where the soft excess component is likely 
to have maximum effect. In the following sections we shall exten- 
sively investigate the effect of a SEDs, given by Equations|4]and[5] 
on the nature of the WA, as a function of the shape of the accretion 
disk component parametrised by the value of Tn, and the strength 
of the soft excess feature parametrised by rj or Fc and 77'. 



3 STABILITY CURVES ANALYSIS 

Studies of WA variability in response to continuum changes 
show that it is reasonable to assume the WA to be in ion- 
ization and thermal e quilibrium as observed for NGC 98 5 
jKrongold et akl l2005al). NGC 3783 jKrongold et al. I l2005bl) . 
NGC 5548 jAndrade- Velazquez et al. I2OI0I) . The thermal and ion- 
ization equilibrium is governed by heating due to photoionization 
and cooling due to line emission and collisional recombination (ra- 
diative and dielectronic). 

We model the WA as an optic ally thin, plane parallel slab 
of So lar metallicity g as (as given by Allende Prieto et al. 1 l l200ll 
2002h for C and O bvlHolwegej ( l200lh for N. Ne. Mg. Si and by 
Grevesse & Sauval ( Il998h for the remainder of the first thirty ele- 
ments) with column density A'^h ~ W^^ cm~^. The absorbing gas 
is assumed to be illuminated by an ionizing continuum given by 
Equation|4]or Equation |5] and the ionization state of the gas can be 
described by specifying the ratio of the ionizing photon flux to the 
gas density through an ionization parameter (see Equation[6j. 
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Figure 2. Stability curves demonstrating tlie advantage of using as compared to The dotted-and-daslied stability curves in both the panels are drawn 
using the broken power-law continuum (Equation[2), whereas the soHd ones are generated for a SED given by Equation|4]with r) = 0.0 and Ti„ = 10 oV. 
Curves drawn with ^ (left panel) are separated in the phase space even if the WA in both the cases have similar nature. On the other hand, the curves drawn 
using the alternative definition, ^p, of the ionization parameter (right panel), are almost identical showing the true nature of the WAs. See text in Section [3T1I 
for further details. 



3.1 Normalising the ionisation parameter 

The photoionization state of the WA can be parametrised by the 
ionization parameter which is the ratio of the ionizing photon flux 
to the gas density (Tarter, Tucker & Salpeter, 1969); 



du [erg cms "'^],(6) 



hi/0 = 13.6 cV 



where Lion is the luminosity between 13.6 eV — 13.6 keV (i.e 1 — 
10'' Rydberg). Hence ^/T ~ L/pR^, p being the gas pressure in 
the WA. 

This definition has a certain drawback which is demonstrated 
by the stability curves on the left panel of Figure |2] The dotted- 
and-dashed curve is for a gas illuminated with the broken power- 
law ionizing continuum (Equation |2j and the solid curve is gen- 
erated with a more realistic SED given by Equation |4] having an 
accretion disk component with Tin = 10 eV, but no soft excess 
(r] = 0.0). From Figure [T] we can see that the SED including the 
accretion disk component (the solid black line) has ~ 50 times 
more photons at 13.6 eV (1 Rydberg) than the dotted-and-dashed 
spectra (Equation |2j, although they have similar flux at 100 eV. 
This would result in very different values of ^ in the two cases 
but very similar WAs, since the WA properties are determined by 
the photon distribution in soft X-ray {E > 100 eV) and not by 
photons with energy E « 100 eV. The kinks in any stability 
curve are a result of the interplay between the various heating and 
cooling agents resp onsible for main taining the gas at a state of 
thermal equilibrium jchakravortv et al. 20 08. 20Q9. ). The stability 
curves on the left panel of Figure|2]have almost identical shape for 
4.4 < logT < 6.5, indicating that they have WAs with identical 
thermal properties and state of ionization. However, they are sepa- 



rated from each other, by a mere horizontal shift in Iog(5/T) be- 
cause of the different ^ values predicted by the two different SEDs. 
Thus, ^ is a relatively poor parametrization of the nature of the 
WA. Such a problem with the standard definition of the ionization 
param eter has been ac knowledged by other authors as well. For ex- 
ample, IChelouche & Netzei, (200^ and references therein) use an 
ionization parameter Ux which considers the ionizing flux only be- 
tween 540 eV to 10 keV. 

To circumvent this problem we use a normalisation scheme, 
appropriate for this paper, as described in the following. All the 
SEDs considered in this paper (except for ones in Section lSTt have 
a common power-law in addition to the diskblackbody and/or the 
soft excess component. Hence we define 
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where Lp is a constant luminosity due to the soft X-ray power- 
law component f{u) ^ v~°' (with a = 0.8) in the energy range 
E = hv = 13.6 eV - 200 keV. Thus, ^p is the ratio of the 
ionizing flux, due only to the power-law component in the energy 
range 13.6 eV to 200 keV, to the gas density. For any SED given 
by Equations |2] |4] or [5] there would be a unique factor ^/^p for 
the same riH R? ■ The stability curves corresponding to the ionizing 
continua given by Equations (2] |4] or |5] are calculated by CLOUDY 
taking the entire SEDs into account, following which their x-axis 
(log(^/r)) is divided by the unique (to each SED) factor ^/^p gen- 
erating a normalised stability curve in the log T - log(^p/r) plane. 
The advantage of such a normalisation is shown in the right panel of 
Figure |2] where the stability curves overlap closely, demonstrating 
the true physical scenario of two very similar WAs. Note that these 
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are the same stability curves which are seperated by the horizon- 
tal shift in the left panel, as discussed above. This normalisation 
scheme is used for all subsequent figures in the paper except for 
Figure [9] In addition to the SEDs which determine the physics of 
the respective stability curves, the top panels of Figures [3l4l and l5l 
also include the constant power-law component (magenta, long- 
and-short-dashed curve) which is used to normalise the stability 
curves so that they can be plotted in terms of log (^p/T). 

It is to be noted that while calculating the thermal and ioniza- 
tion properties of the WA, CLOUDY uses the entire SEDs, including 
all three spectral components (namely the disk blackbody, the soft 
excess and the X-ray power-law), described by Equations |4] and [5] 
The introduction of and normalisation of the stability curves is 
only for the convenience of demonstration of the stability curves 
and their associated properties. Use of ^p merely introduces an 
overall horizontal shift in these distributions. Further, note that any 
given range in log(^/r) would correspond to an exactly equal 
range in log(5p /T) and vice versa. 
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3.2 Disk-blackbody 

The top panel of Figure [3] shows the ionizing continua which in- 
clude the disk blackbody and the power-law components (77 = 0, 
Equation 0. As Tin is increased from 10 eV to 30 eV, the peak 
of the flux distribution moves from ~ 30 — 90 eV. Photons at 
i? < 90 eV affect only the lower ionization species of the absorb- 
ing gas. The resulting WA should show differences only in the na- 
ture of the lower ionization phases, and should have very similar 
higher ionization states. We shall discuss this issue further in Sec- 
tion|4] 

Corresponding stability curves are shown in the bottom panel 
of Figure [3] A hotter disk blackbody lowers the Compton temper- 
ature branch and increases the temperature of the bottom branch. 
However, in the 5 < logT < 6.5 K temperature range the WAs, 
the stability curves are independent of the value of Tin, i. e. , unaf- 
fected by the spectral component due to the accretion disk. How- 
ever, the lower temperature (log T < 5) part of the stability curves 
undergo remarkable changes with the increase in r,„ ; for the same 
value of ^p, the gas attains a higher temperature, the stability curve 
becomes more stable and any possibility of multi-phase existence 
with higher ionization states is lost. See Section |6] for further dis- 
cussion on this issue. 

3.3 Soft excess 

3.3.1 Blackbody soft excess 

The top panel of Figure |4] shows the ionizing continua given by 
Equation |4] with values of 77 = 0.0, 1.0, 3.0, representing increas- 
ing strengths of the soft excess component for a fixed value of the 
blackbody temperature T^e = 150 eV. The corresponding stabil- 
ity curves in the bottom panel show a remarkable enhancement of 
the stable log(^p/r) range from ~ 0.20 through 0.41 to 0.63 
dex for the 10'' K absorber, as the strength of the soft excess is 
increased from ry — 0.0 through 1.0 to 3.0. The 10*^ K phase re- 
mains unchanged through the variation of 77. In the 77 = 3.0 case 
however, the transition from the ~ 10^ to 10® K phase is smooth 
i.e. with no distinct unstable region separating the two temperature 
regimes. Such a result implies that stronger the soft excess compo- 
nent in the spectra, greater is the probability of finding a ~ 10^ K 
absorber, since the stable region then becomes less susceptible to 
luminosity variations in the AGN. However, the possibility of the 
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Figure 3. The SEDs in the top panel are constituted by two spectral com- 
ponents, a X-ray power-law and a EUV disk blackbody. The curves are for 
different values of the temperature Ti„ of the innermost ring of the accre- 
tion disk, a = 0.8 and aox = 1-2 for all the SEDs. We have included 
the constant power-law (magenta, long-and-short-dashed curve) which is 
used to normalise the stability curves so that they can be plotted in terms 
of log(5p /T). Important energy values including the upper energy cut-off 
{Emax = 200 keV) and the range of definition for aox (2500A and 2 
keV) have been marked and labeled. The bottom panel shows the stability 
curves corresponding to the ionizing continua shown in the top panels. 



10^ K gas being in pressure equilibrium with the other WA phases 
is significantly reduced. We shall return to these points in Section[6] 
To find a typical value of r) implied by X-ray observations of 
AGN we refer to two examples, namely a typical Seyfert 1 galaxy 
NGC 5548 and a narrow l i ne Sey fert 1 galaxy IRAS 13349-1-2438. 
lAndrade- Velazquez et al. I ( I2OI0I) analysed the 800 ks Chandra 
grating spectra of NGC 5548 and their best fit parameters for the 
soft X-ray continuum is a power-law with a = 0.6 and a blackbody 
with Tse = 110 eV. The relative normalisation of the two com- 
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Figure 4. Top panel : SEDs wliicli liave a soft excess component de- 
scribed by a blackbody distribution (see Equation |4) with temperature 
Tse = 150 eV. The different SEDs are for different values of r], i. e. , for 
different strengths of the soft excess component, a = 0.8 and apx = 1-2 
for all the SEDs. The magenta, long-and-short-dashed curve is the same as 
in the top panel of Figure [5] We have included the energy values 1300A 
and 120 eV to guide the eye in the context of observations of Mrk 478 (see 
text in Section D.3.2) Bottom panel : The stability curves corresponding to 
the ionizing continua shown in the top panels. With the growing strength of 
the soft excess component, there is significant increase in the range of 
over which thermally stable warm gas at ~ 10^ K can exist. 



ponents results in 77 ~ 1.2. Similarly. iHolczer et al. 1 toot have 
analysed the 300 ks Chandra data for IRAS13349+2438 and have 
reported the spectral parameters for the best-fit ionizing continuum 
comprising of a X-ray power-law with q = 0.9 and a blackbody 
with Tse = 105 eV. From their results we calculate 77 to be ~ 2.58. 
In both cases the values of i] is well within the range considered 
here (0 - 3). 

The Tse, observed for a large number of type I AGN, lies in 



the small range 100 — 200 eV. We have checked whether the influ- 
ence of the soft excess component is modified when the blackbody 
temperature is varied across this range. Results showed that the sta- 
bility curves remain unaffected by such variations, e.g. the range in 
log(^p /T) for the stable 10''' K gas changes by merely 0.06 dex as 
Tse varies from 100 eV to 200 eV. 



3.3.2 Comptonised soft excess 

The blackbody soft excess sometimes turns out to be inadequate 
in explaining the observations. For example, the SED for Mrk 
478 may be 1-3 times brighter in vLu at E ^ 120 eV than it 
is at 1300 A (~ 10 eV) l lGondhalekar et al. II 19941 : [Marshall et all 
fl99^ . We find that a blackbody soft excess cannot satisfy such a 
condition even with a relatively high normalisation of 7; = 3.0 
(top panel, Figure|4]l. We want to investigate if the alternative de- 
scription of the soft excess component, namely the soft excess due 
to thermal Comptonisation, can account for a SED similar to that 
seen in Mrk 478, and to see what are its effect on the WA. 

Spectra generated by Equation |5] with fixed values of Tc = 
150 eV, Tin = 10 eV and 77' — 1.0 are drawn in the top 
panel of Figure |5] for different values of Te- For comparison, 
we have also drawn the continuum with no soft excess (Figure |5] 
solid curves). The continuum with Fc — 2.0 satisfies the con- 
ditions required by o bservations of Mrk 478 as mentioned by 
iKorista. Ferland & Baldwi n (1997). 

The influence of a thermal Comptonisation soft excess com- 
ponent on the stability curves is shown in the bottom panel of Fig- 
ure |5| We see again, that the range of log(^p/r) for the stable 
lO"" K branch is increased in width from ~ 0.2 dex to 0.48 dex 
as 77' goes from to 1.0, independent of the value of F, facilitat- 
ing the presence of a stable WA at these temperatures and that the 
different phases of the WA lose the possibility of existing in pres- 
sure equilibrium as the strength of the soft excess component is 
increased. 

The thermal Comptonised and the blackbody models of the 
soft excess result in similar WAs. For example, 77 = 1.0 curve 
(dashed; lower panel of Figure|4|( have stable ranges of log(^p/r) 
at 10^ and lO" K similar to that of the 77' = 1.0 curves (the dot- 
ted and the long-dashed for F = 1.7 and 2.0 respectively; lower 
panel of Figure|5](- The stability curves with 77' — 1.0 do have a 
slightly smoother rise from the ~ 10^'^ to 10* K with no distinct 
intermediate unstable phase, because they have more flux in the 
i5 = 10 — 100 eV range for the same soft excess strength. 



4 CAUSE OF EXTENDED STABLE REGION : HEATING 
AGENTS AND ION FRACTIONS 

We have seen that the disk blackbody and the soft excess radia- 
tion at < 0.5 keV significantly influence the stability curve, im- 
plying that these components are important for shaping the ionic 
and thermal state of the WA (SectionjS}- Changes in the disk black- 
body affect the low temperature arm (logT ~ 4.5) of the stabil- 
ity curve, whereas the soft excess influences the absorbing gas at 
log T ^ 5.0 K. In this section we examine the important heating 
agents and the distribution of ion fractions of the signature ions in 
the WA along the stability curve. All the ionizing continua con- 
sidered in this section are given by Equation |4| The effect of the 
thermal Comptonised soft excess gives the same qualitative results 
as the blackbody soft excess modeled by a blackbody (Equation |4| 
see ). 
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Figure 5. Top panel : Ionizing condnua witli soft excess component due to 
tliermal Comptonisation given by Equation \5\ witli fixed values of Tc = 
150 eV, Tin = 10 eV and rj' = 1.0. The different curves are for different 
values of Fc, namely 1.7 and 2.0. For comparison, we have also drawn the 
continuum with no soft excess (solid curve). With a = 0.8 and apx = 
1.2, the accretion disk component and the X-ray power-law component are 
same for all the SEDs. The magenta, long-and-shoi1-dashed curve is the 
same as in the top panels of Figures[3]and|4] We have included the energy 
values 1300A and 120 cV to guide the eye in the context of observations 
of Mrk 478 (see text in Section 13.3.2) Bottom panel : The stability curves 
con'esponding to the ionizing continua shown in the top panels. 



4.1 Heating 

CLOUDY works by dividing a ga.s into a set of thin concentric shells, 
referred to as 'zones' which have thicknesses that are small enough 
for the physical conditions across them to be nearly constant, main- 
tained by continuously adjusting the physical thicknesses of these 
shells. For this sub-section, only, we have constrained the code to 
perform a single zone calculation for numerical convenience as we 



are not concerned with quantitative rigor, but want to understand 
the qualitative trends. 

Figure |6] shows the influence of these two components on the 
heating fraction Aff of H+", He+\ Fe and O as a function of 
(,p/T. All other ions and/or elements contribute less than 10% to 
AH. 

On the left the solid and the dotted-and-dashed curves re- 
spectively correspond to a WA ionized by an SED with Tin = 
10 and 30 eV (see Eqn|4l. In both cases 77 = so that there is 
no soft excess component in these SEDs, so any changes in the 
heating fractions are entirely due to the variation in the accretion 
disk spectrum. 

On the right panels the solid and the dotted-and-dashed curves 
correspond to 7; = and 3 respectively, while Tin = 10 eV. 
Hence changes in the heating fractions in the right panels are due 
only to the varying strength of the soft excess component. 

The top panels (labeled A) of Figure [6] show the region of the 
stability curves where the variation of these spectral components 
has maximum effect. 

4.1.1 Hydrogen and Helium 

Panels B and C of Figure [6] show the distribution of AH respec- 
tively for neutral Hydrogen H"*"" and singly ionised Helium He"^^. 
In the left panels we see that there is no significant difference in 
the distributions as Tn increases from 10 to 30 eV. For both H+" 
and He+^, AH is considerably less for r) = 3.0 than for 77 = as 
seen in the right panels. Thus in either case (of the disk blackbody 
or the soft exces.s) H^" and He"'"^ can be ruled out as the cause of 
extra heating causing the increase in the temperature of the stability 
curves in the range of log(^p /T). 

4.1.2 Iron 

We find that in the disk blackbody and the soft excess cases the sig- 
nificant heating agents are the various species of iron. It is to be 
noted that the range of AH for Fe (panels DI, D2 & D3) is double 
that for O (panels El, E2 & E3). The AH behaviour of the iron 
ions are affected, both, when we increase Tin from 10 to 30 eV, 
and when we increase rj from to 1. Fe"'"^,Fe^*, Fe"*"', Fe^^" 
and the Unresolved Transition Array (UTA, see below) are the most 
significant heating agents; and their contributions have been added 
to give the heating fraction due to iron (labeled 'Iron') in panels Dl 
through D3. On the left panels D2 and D3 the heating fractions of 
the individual ions of iron and the UTA are shown, respectively, for 
Tin = 10 and 30 eV. Similarly, the right panels D2 and D3 show 
AH for the individual iron ions and the UTA, res pectively, for 
ri = and 3.0. A WA which is overabundant in iron jFields et al. I 
I2005ll2007h and/or illuminated by a SED having photons facilitat- 
ing absorption by iron will be warmer at relatively lower values of 
^p. These two physical scenarios influence different ions of iron, 
which may help to distinguish between the otherwise degenerate 
effects. 

For -4.0 < log(5p/r) < -3.4, where AH for iron is 
most affected by Tn, Fe"*"* and Fe"*"® are the dominant ions (panel 
D2 and D3 on the left of Figure |6ll. The ionization potentials for 
Fe+^Fe+^ Fe+^ and Fe+^° are 151.1, 235, 262.1 and 290.4 
eV respectively. As the accretion disk becomes hotter, with Tn in- 
creasing from 10 to 30 eV, the number of photons with E > 100 eV 
increases by a factor of 27. It is this enhancement in the number 
of photons that influences these ions of iron differently in the two 
cases and hence changes the shape of the stability curves. 
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Figure 6. The principal heating agents which are being influenced by the soft spectral components disk blackbody (panels on the left) and soft excess 
(panels on the right) present in the SEDs of typical AGN. Top panels (A) : The stable part of the thermal equilibrium curves for an optically thin Solar 
metallicity gas illuminated by ionizing continua given by Equation |4] On the left panel the solid and the dotted-and-dashed lines correspond to SEDs having 
Tin = 10 and 30 cV respectively with = in both cases. The solid line in the right panel represents a SED with r] = and the dotted-and-dashed curve 
is for r] = 1.0, where Ti„ = 10 cV for both the curves. We have zoomed in on the regions of the stability curves where they differ most from one to the 
other. The line schemes described for the left and right panels are maintained same for the lower panels in the figure. Lower panels : Fraction ( AH) of the 
total heating caused by the various significant elements and ions. We have demonstrated the effects of all those elements which contribute > 10% to the total 
heating in the same log(^p/T) range as shown in the top panels. Panels B, C, Dl and El show the results for H+'\ He+^, iron and oxygen respectively. 
The 'break-up' in heating fraction of iron are shown in panels D2 and D3 and the same for oxygen are shown in panels E2 and E3. Note that the y-range for 
the panels El through E3 for Oxygen are different from all other panels showing AH distributions because Oxygen never contributes more than 20% to the 
heating unlike other agents, some of which are responsible for even up to 50% of the heating for some values of log(5p /T). Note that the discontinuity in the 
AH distribution of oxygen in the = 3.0, Tin = 10 eV case (dotted-and-dashed curve in panels El and E3) is because the contribution drops below 5% for 
-4.2 < log(Cp/T) < 4.0. 
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The UTA is a blend of the numerous absorption lines aris- 
ing from the iron M-shell ions, Fe'*''' - Fe^'^^, due to their n = 2-3 
(mainly 2p-3d, n being the principal quantum number of the active 
electron) transitio ns which are found to o ccur betw een 16 and 17 
A (730 - 776 eV) jSako et al 11200 iklNetzer 2004 : Krongold et al. I 
l2005al lbl: I Holczer et al. II2007I . and reference therein). The energy 
ranges important for the UTA are likely to be influenced by the soft 
excess component of the SED. This effect is clearly seen by com- 
paring the Aiif distributions of UTA in the right panels D2 and D3 
of Figure [6l AH for iron is far more dominated by the UTA AH 
contribution for t) = 3.0 case in the range —4.4 < log(^p/T) < 
3.6. 

Note that in panels D2 (both left and right), the solid curve 
for iron is multivalued at -3.43 < \og{£,p/T) < -3.42, just 
where the stability curve is multivalued (see panels A), giving mul- 
tiple phases in pressure equilibrium. On the other hand, for the 
Tin = 30 eV and ther; = 3 curves (the dotted-and-dashed sta- 
bility curves in panels A and the corresponding AH distribution of 
iron in panels D3, both left and right) there is no multivalued be- 
haviour for log(^p /T) ~ —3.4. Thus the behaviour of the stability 
curve seems to be driven by the heating due to the different ions of 
iron. 



4.1.3 Oxygen 

Panels El through E3 show the results for the heating fractions 
contributed by Oxygen. Note the factor 2 smaller, y-range for the 
panels El through E3 because Oxygen never contributes more than 
20% to the heating of the gas. 

The left panels show that AH due to Oxygen is lower for 
Tin = 30 eV than that for Tn = 10 eV in the range log(C/r) < 
—3.6. In the same range Of log(^/T), the stability curve for Tn = 
30 eV has higher temperature than the Tn = 10 eV curve (panel A 
on the left). Thus here Oxygen is not acting as one of the required 
heating agents resulting in a warmer absorber. 

However, in the 77 = 3.0 case (right panels) Oxygen becomes 
a significant excess heating agent for log(^p/T) > —3.6 where 
the contributing ions are 0+*^ (OVII) and 0+'^ (OVIII). 

As for iron, the AH distribution due to oxygen for the Tn — 
10 eV, 77 = case, are also multivalued at the same values, namely 
the narrow range —3.43 < log(^p/r) < 3.42 (panels E2, both 
left and right), although not as strongly as that for iron (panels D2, 
both left and right). Moreover, at log(^p /T) --^ —3.4, it is iron and 
not oxygen which is the dominant heating agent. However, oxygen 
is the dominant heating agent at higher values log(^p/T) ~ —3.2, 
and has multivalued AH distribution which drives the multi-phase 
nature of the solid stability curves at those values of log(^p /T). 

The AH distribution of oxygen for the Ti„ = 30 eV, 77 = is 
shown in the left E3 panel (dotted-and-dashed curve). For — 3.8 < 
log(^p/T) < —3.3 the oxygen AH is about 2 to 4 times lower 
than the iron AH (left panel D3). The oxygen AH dominates only 
beyond log(^p/r) > —3.3 where it is multivalued and is respon- 
sible for the multivalued nature of the stability curve (left panel A) 
at the same values of log(^p/r). 

For the 77 = 3.0, Tn = 30 eV case (right panel E3) as well, 
the oxygen AH contribution is lower than the iron AH (right panel 
D3) until log(^p/r) > —3.6 beyond which the oxygen AH dis- 
tribution determines the behaviour of the stability curve. Note that 
the stability curve (dotted-and-dashed in right panel A) is not mul- 
tivalued at log(5p/r) > —3.6 because the AH for oxygen is not 
multivalued. 

Thus, at the higher ionization (log(^p/r) ~ —3.2) phase of 



the WA, the nature of the stability curve seems to be driven by the 
heating due to the oxygen ions, predominantly O^^ (OVIII). 



4.1.4 Zero oxygen and iron abundance 

We have further illustrated the importance of iron and oxygen as the 
key players in determining the thermal properties of the WA using 
Figure[7] The distinction between the left and the right panels are as 
that in Figure|6j the upper panels being the same stability curves as 
in Figure[6]panels A, but for slightly different ranges of log(^p /T) 
and wider ranges of log T. 

Comparing the stability curves for Solar metallicity gas (top 
panels), which have pronounced kinks in the temperature range 
4.2 < logT < 6.5, with the relatively featureless ones for gas 
with zero O and Fe (bottom panels), we see that the detailed shape 
of the curves in the temperature range ~ 4.2 < log T < 6.5 are 
go verned by the atomic inte ractions due to oxygen and iron (also 
see IChakravortv et al. II2009I) . As a result, the SED induced differ- 
ences between the solid and the dotted-and-dashed stability curves 
in the top panels are significantly smoothed out in the bottom pan- 
els. 

For —3.7 < log(^p/r) < —3.3 the rise in temperature of 
the Solar metallicity stability curves in the top left panel shows 
a temperature difference of 0.4 dex from one curve to the other. 
This difference is reduced to 0.2 dex when iron and oxygen are not 
present in the gas (lower left panel) suggesting that a hotter accre- 
tion disk (Tin ~ 30 eV) mainly affects ions of iron and oxygen 
which results in a hotter WA. 

The temperatures of the Solar metallicity stability curves in 
the top right panels are remarkably different, by '-^ 1 dex, for 
—4.2 < Iog(^p/r) < —3.3. These differences are significantly 
reduced to 0.3 - 0.5 dex when iron and oxygen are removed from 
the gas (lower right panel). This indicates that atomic interactions 
due to these elements are affected by the strength of the soft ex- 
cess component. Unlike the accretion disk case, however, the sta- 
bility curves in the lower right panel retain some of their differences 
in WA temperatures suggesting that other elements have a role as 
heating agents in the soft excess case. The detailed investigation of 
this effect is beyond the scope of this paper, but will be attempted 
soon in our future publications. 



4.2 Ion fraction 

Absorption lines and edges of OVII, with ionisation potential (IP) 
of 0.74 keV and OVIII, with IP of 0.87 keV, are often prominent 
signatures of the WA in soft X-ray spectra. Hence the column den- 
sities of these two ions are often considered as important observ- 
able parameters for WA states (see introduction for references). At 
higher values o f ^p, absorption by various ions of silicon become 
important, e.g. iNetzer et al.1 j2003h show that SiXIV (IP of 2.67 
keV) has a significant column density conesponding to the high 
temperature component of the WA. We choose these three ions to 
study the variation of their ion fraction as a function of ^p as the 
spectral components change in the ionizing continuum. 

The ion fraction of the ion is the fraction of the 

total number of atoms of the element X which are in the i"^ state 
of ionization. Thus, 



I{X+') = 



N{X+') 



fix) Nu ' 

where N{X'^^) is the column density of the X^^ ion and f{X) 
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Figure 7. Top Panels : The compaiisoii between the stabiHty curves for Solar metalhcity gas, when ionized by SEDs having accretion disks with different 
temperatures Tin = 10 and 30 eV (left) and when the ionizing continua have different strengths of the soft excess component rj = and 3.0 (right). Bottom 
panels : The same ionizing continua are used to illuminate a gas for which iron and oxygen have been removed. The relative differences in the stabiUty curves 
are almost entirely removed in the lower left panel, and are significantly reduced in the lower right panel, indicating that these spectral components strongly 
influence the ions of iron and oxygen. 



n{X)/nu is the ratio of the number density of the element X to 
that of hydrogen. 

The top panels of Figure[8]show the distribution of log T with 
respect to log^p for the stable WA phases. The gaps in the lines 
correspond to the range of ionization parameter over which the gas 
is thermally unstable. The lower panels show I(OVII), I(OVIII) and 
I(SiXIV). 

As before, left panels in Figure [8] correspond to the changes 
in the diskblackbody. For Tin ~ 10 eV we can see three distinct 
regions of separated by intermediate ranges of thermally unsta- 
ble solutions. These regions are respectively dominated by OVII, 
OVIII and SiXIV. However, when Tin is raised to 30 eV, the sta- 
bility curve becomes stable all the way up to log T = 5.3 and the 
OVII dominated phase gradually merges in to the OVIII dominated 
gas, as demonstrated by the uninterrupted dotted-and dashed line in 
the top left panel. Thus there seems to be a continuous distribution 
of possible thermal solutions resulting from the greater production 
of OVII (second panel from top) which, in its turn, is a result of the 
enhanced number of i5 > 100 eV photons coming from a hotter 



accretion disk with Tn ~ 30 eV. The change in the shape of the 
disk blackbody does not change the ion fraction distribution of the 
high ionization species like OVIII or SiXIV (the two lowest pan- 
els) or the stability curve for logT > 5 and log^p > 1.7. Thus 
the higher temperature (log T > 5) phase of the ionized absorber 
remains unaffected. 



The effect of the soft excess component is demonstrated in the 
right panels of Figure [8] In this case the increase in the strength 
of the soft excess component from = to 77 = 3.0 makes the 
absorber thermally stable all the way up to temperatures of log T ~ 
6. The third panel from the top shows that increase in the strength 
of the soft excess component mainly influence the occurrence of 
OVIII. The fraction of OVII is reduced for rj — 3.0 because of the 
facilitated production of OVIII at lower values of log^p. This is 
caused by increase in the number of ^ 100 eV — 2 keV photons by 
a factor of ~ 7.5. The ion fraction distribution of SiXIV however, 
remains unaffected (bottom panel). 
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Disk blackbody Soft Excess 




Figure 8. The ion fractions of some of the relevant ions as a function of ionization parameter; the left panels show how this distribution is influenced by the 
disk blackbody spectral component and the right panels show the same for the soft excess. Top panels : The distribution of temperature as a function of is 
plotted for the thermally stable WA. The gaps in the curves correspond to the range of over which the stability curve is thermally unstable. An increase in 
accretion disk temperature Tin from 10 to 30 cV makes the gas stable at log T ~ 5. When the strength of the soft excess is increased from r; = to 3.0, all 
the phases of the WA become thermodynamically stable up to a temperature of log T 6.3 Lower panels : Ion fraction for various ions (as labeled) important 
for the WA states. The line-style scheme followed here is the same as in Figure|6]and is also labeled in the second panel from the top. See text for the physical 
impUcations of the ion fraction distribution. 



5 OTHER SIGNIFICANT PHYSICAL PARAMETERS 
FROM EARLIER STUDIES 

IChakravortv et al. had studied the properties of the WA as a 

function of the X-ray spectral index a and the chemical abundance 
of the gas. In this subsection we summarise some of those results 
which are relevant in the context of the soft spectral components of 
typical AGN continua and their effect on the WA. 



5.1 X-ray spectral index a 

IChakravortv et al. I ( l2009l Section 5.1 and Figure 3) showed that 
there was no thermally stable absorber for logT > 4.4 for a flat 
spectral index of a = 0.2 (in Equation[2l aox = 1.2). Figure|9] 
shows the stability curve using Equation [2] with a — 0.2 (dotted, 
magenta line), which corresponds to an ionizing continuum with 
no soft excess component, shows no appreciable stable WA phase 
at T ~ 10^ K. In Section g] of this paper we have seen that the 
presence of soft excess component in the AGN spectra facilitates 
the existence of thermally stable gas at logT ~ 5.0. However, 
throughout we have used a = 0.8. It is thus crucial to check if 
the soft excess component can result in stable 10^ K gas even if 
a ~ 0.2. 



The solid, the dotted-and-dashed and the dashed curves in Fig- 
ure [9] show the comparison between the stability properties of the 
WA ionized by SEDs given by Equation |4] with Ta^ ~ 150 eV, 
ri — 3.0 and r,„ — 10 eV but various values of a (0.8, 0.5 and 
0.2). The range, in log(^/r), of stable WA phase at log T ~ 5 de- 
creases as a function of decreasing value of a. However, even for 
Q as low as 0.2 we still see a stable phase at log T 5. Thus the 
appearance of a stable state at logT ~ 5.0 for the a = 0.2 curve 
can be attributed to the presence of a strong soft excess component 
present in the ionizing continuum. However, if a < 0.2, even with 
contributions from a strong soft excess component, there is no ev- 
idence for stable gas at any higher temperature than logT > 5.3 
(refer to the dashed, blue curve). 

Notice that in Figure |9] the stability curves are drawn using ^ 
instead of ^p. The normalisation scheme for the ionisation param- 
eter discussed in Section [3T| is not appropriate while comparing 
SEDs with different values of the X-ray spectral index a. 

5.2 Chemical abundance of the absorber 

We discuss the role of the elemental abundances of the absorbing 
medium in IChakravortv et al. I ilOO^ . Here, we revisit the super- 
Solar metallicity results, but, with more realistic AGN continua 
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Figure 9. Thermal equilibrium curves generated using SEDs given by 
Equation |4] with a strong soft excess (Tse = 150 eV, rj = 3.0 and 
Tin = 10 eV) for a range of X-ray slope a = 0.8 0.5 and 0.2 (respec- 
tively drawn with the solid, the dotted-and-dashed and the dashed lines). 
The stable phase at logT ~ 5.0 the range of log(^/T) progressively de- 
creases as a decreases from 0.8, to 0.5 and then to 0.2. For comparison we 
have used Equation|2]with a = 0.2, which corresponds to an ionizing con- 
tinuum which has no soft excess component, and drawn the dotted stability 
curve which is consistent with no stable WA phase at log T ~ 5.0. 



given by Equation|4](a = 0.8, rj = 3.0 and Tin = 10 eV). The sta- 
bility curves are shown in Figure[TO]for WA abundances from Solar 
( Zp)) to 9 Zfl in steps of 2. The qualitative results are the same as 
in lChakravortv et al. I lioosl) . 

On the low temperature (logT < 4.5, corresponding to OVII 
like ions) arm of the stability curves, super-Solar metallicity results 
in a cooler absorber for the same (or 5) values, with an increase 
in the range of over which we get stable WA. Thus super-Solar 
metallicity opposes the role of increased Tin (see Figure[3]in Sec- 
tion[T2ll. 

For the intermediate temperature arm of the stability curve 
(logT ~ 5.0, corresponding to OVIII like ions) super-Solar abun- 
dance decreases the range of stable WA, thus opposing the influ- 
ence of the increase in rj in the ionizing continuum (see Figure|4]in 
SectionlU. 

On the highest temperature arm of the WA (log T ~ 6.0, cor- 
responding to SiXIV like ions) super-Solar gas tends to be hot- 
ter, and is thermodynamically stable for a larger range of values of 
log(5p /r). Thus the fine tuning of values to detect such species 
of ions is relaxed if the absorber has super-Solar abundance, and the 
probability of detecting these high ionization species is increased. 
In this regime, the influence of high metallicity is degenerate with 
the influence of increasing a in the ionizing continuum, but only if 
a > 0.8 (see Figure[9]in Section lSTb . 



6 MULTI-PHASES WARM ABSORBER 

The stability curves we have derived in the previous sections of- 
ten contain segments that allow phases at different temperatures to 
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Figure 10. Stability curves for absorbers with Solar and super-Solar abun- 
dance from Zq to 9 Zq in steps of 2. The ionizing continuum has a strong 
soft excess (Equation|4]with ; a = 0.8, rj = 3.0 and Tin = 10 eV). 



occur at similar pressure £,p/T indicating the possibility of two 
or more distinct WA phases within the same medium to exist at 
near pressure equilibrium. However, variations in some of the phys- 
ical parameters can also result in the possibility of multi-phase ab- 
sorbers being lost. 

Observationally, distinct phases in pre ssure equilibrium have 
been derived in a number of cases of WAs dKroneold et al. II2003I; 
Netzer et al. ' '2003'; 'Chelouche & Netzed l2005l : iKrongold et al. I 



2007,; Andrade- Velazquez et al. 20_10|), while in other cases it has 
been cla imed that there is a continuous range of ionization pa- 
rameter l lOgle et al. ll2004l ; ISteenbrugge et al. ll2005l) . Since a wide 
range of physical parame ters have been studied in this paper and 
IChakravortv et al. I ( l2009t) for their effects on the WA, we can now 
form a clearer picture of which physical conditions favour the 
multi-phase absorbers. 

The disk blackbody component from the accretion disk affects 
the lower temperature part of the stability curve (Sections[3]and|4]l. 
Figure [3] shows that the stability curve (solid line) with Tn ~ 
lOeVhas a slim range —3.34 < log(^p/T) < —3.31 over which 
gas at low temperatures log T < 4.3 can be in pressure equilibrium 
with the absorber at intermediate temperatures logT < 5.0. How- 
ever, as Tin is increased, all possibilities of such multi-phase struc- 
ture is lost, instead the curve becomes stable in the intermediate 
temperature range supporting a continuous distribution of allowed 
pressure and temperature in this phase space. Thus cooler accretion 
disks support a multi-phase WA, although for a narrow range of 
pressure. 

An ionizing continuum having no soft excess (solid line in Fig- 
ures |4] and [5]l allows for a multi-phase between logT ~ 5 and 6. 
However, with the increase in the strength of the soft excess com- 
ponent the range of log{^p/T) over which the possibility of multi- 
phase exists reduces from 0.09 dex for 77 = to 0.006 dex for 
= 1.0 and to for?? = 3.0. 

IChakravortv et al. I ( l2009h had emphasised that super-Solar 
abundances increase the range of log(5/T) with the possibility 
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of multi-phase WA (also see iKomossa & Mathuj[200ll) . For ex- 
ample, as the metallicity of the gas was increased from 1.0 to 3.0 
to 5.0, A45[log(C/r)] increased from 0.012 to 0.08 to 0.1, where 
A45[log(^/T)] gives the range of common values of log(^/r) 
over which the WA has stable phases at both ~ 10^ and ~ 10'' K. 
Beyond a metallicity of 5Zq, A45 [log(^/r)] was found to drop 
because the range of log(C/T) for the ~ lO'^ WA was found to 
decrease with the increase in metallicity and became zero for a gas 
with metallicity 9Zq . 

Similar trends are retained for the continua in this paper. As 
the metallicity of the gas is increased, keeping the ionizing con- 
tinuum the same, the stable phases increasingly align along the 
same values of log(^p/r) (Figure [TOt. strongly suggesting the 
presence of multi-phase WA. A45[log(5p/T)] increases from to 
0.02 to 0.04 to 0.045 as metallicity is increased from Solar to 3, 5 
and 7 times Solar. Beyond this high metallicity, A45[log(^p/r)] 
becomes constant. There is a even more significant increase in 
A56[log(^p/T)] from to 0.03 to 0.08 to 0.11 to 0.13 with the 
increase in metallicity from Solar to 3 , 5,7 a nd 9 times Solar. 

For NGC 3783, iNetzer et al.] ( l2003h found three distinct 
phases of the WA in near pressure equilibrium. In Figure [TO] 
we see that only for Z = 9Z0 there is a very narrow range 
of 0.03 dex in which all three phases of the WA, at tempera- 
tures ~ 10'',10''^ and lO'^ K, are stable in the range -3.69 < 
log(^p/T) < —3.66. Interestingly, the stability curve for a much 
lower metallicity, Z = 5Z0, gas but ionized by a continuum with 
weaker soft excess (rj = 1.0, not shown in Figure (Toj shows the 
three phases in pressure equilibrium for a larger range of 0.08 dex 
for —3.48 < log(^p/T) < —3.4. Thus, it is easier (relatively 
lower super-Solar metallicity) to have 3 phases in pressure equi- 
librium if the strength of the soft excess is low er. In all the physi- 
cal co nditions of the WA, that were studied bv lChakravortv et al. I 
( l2009h . three phases of the WA were found to coexist in pres- 
sure equilibrium for still lower met allicity > SZp, which br ings 
it down to observed values (e.g. see [Fields et al. |[2005l . l2007l for 
MRK 279). It is to be no ted, however, that the continua used by 
IChakravortv et al. I ilOO'^ had no soft excess component in them. 



7 CONCLUSION 

We have examined the effect of spectral energy distributions in- 
cluding hot disks and soft excess spectral components in the energy 
'blind spot' between 13 - 100 eV (Equations |4] and [Sj on warm ab- 
sorber. We investigated whether the spectral shapes can, in turn, be 
constrained from the observed properties of the warm absorber. 
We summarise our results as follows : 

• The maximum temperature of the accretion disk component 
(see Equation [4] Section [2!2t strongly affects on the low tempera- 
ture (logT < 4.5) arm of the stability curve, the thermal proper- 
ties of which are largely decided by ions Fe+^ to Fe+^'' and 0+^^ 
(OVII). This phase of the WA becomes hotter and thermodynami- 
cally stable over a larger range of ^ when the ionising continuum is 
hotter (Tin = 30 eV) compared to Ti„ = 10 eV. Thus hotter accre- 
tion disks are more likely to produce WA phases characterised by 
ions having similar ionisation potentials to OVII. However, the pos- 
sibility of this phase of the WA being in pressure equilibrium with 
the higher temperature phases is eliminated in the hot disk cases 
as A45[log(C/r)] decreases from 0.03 for Ti„ = 10 eV to for 
Tin = 20 and 30 eV. Changes in the accretion disk temperature 
however, do not affect the higher temperature arms of the stability 
curve. 



• The thermal properties of the intermediate temperature 
(logT ~ 5.0) branch of the stability curve are essentially deter- 
mined by atomic interactions due to the different ions of iron and 
OVIII (IP = 0.87 keV). We find that with the increase in the relative 
strength rj of the soft excess component this 10"" K stable phase of 
the WA spans a much larger range of log(^/r) (by 0.4 dex from 
rj — to 77 — 3.0) thus increasing the probability of finding WA 
phases characterised by 'OVIII like' ions. 

• For a gas whose chemical composition is devoid of iron and 
oxygen, changes in the accretion disk spectral component bring 
about no change in the properties of the WA, and changes due to 
the variation in the strength of the soft excess component are sig- 
nificantly reduced. 

• The highest temperature arm (logT ~ 6.0) of the stability 
curve which is characterised by the ions having similar IP to that 
of SiXIV (IP = 2.67 keV), is left almost unaffected by any changes 
either in the accretion disk component or in the soft excess compo- 
nent. 

• An AGN continuum with a flat X-ray slope will not produce 
high IP ions like SiXIV, and may show signatures of lower IP ions 
like OVIII (IP = 0.87 keV) only if a sufficiently strong soft excess 
component is present in the SED. 

• The metallicity of the gas plays an important role in determin- 
ing the multi-phase nature of the warm absorber. The possibility 
of the lO'' and the lO'^ K phases of the warm absorber occurring 
at similar values of log(^/T) is increased if the absorber is super- 
Solar in abundance. Similarly, the chances of having pressure equi- 
librium between the 10^ and the 10* K phases are also increased 
by the super-Solar metallicity of the gas. However, all three phases, 
together, are found to be in pressure equilibrium only if Z > 3Zq, 
and the required abundance to achieve this is increased if the ioniz- 
ing continuum has stronger soft excess. 
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